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Abstract 

It has long been known that complex balanced mass-action systems 
exhibit a restrictive form of behaviour known as locally stable dynam- 
ics. This means that within each compatibility class C Xo — the forward 
invariant space where solutions lies — there is exactly one equilibrium 
concentration and that this concentration is locally asymptotically sta- 
ble. It has also been conjectured that this stability extends globally 
to C Xo . That is to say, all solutions originating in C Xo approach the 
unique positive equilibrium concentration rather than dC Xo or oo. To 
date, however, no general proof of this conjecture has been found. 

In this paper, we approach the problem of global stability for com- 
plex balanced systems through the methodology of dividing the pos- 
itive orthant into regions called strata. This methodology has been 
previously applied to detailed balanced systems — a proper subset of 
complex balanced systems — to show that, within a stratum, trajecto- 
ries are repelled from any face of IR> adjacent to the stratum. Several 
known global stability results for detailed balanced systems are gener- 
alized to complex balanced systems. 
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1 Introduction 



Chemical reaction modeling is rapidly becoming a topic of great interest 
in areas such as systems biology, atmospherics, pharmaceutics, industrial 
chemistry, etc., where mathematical tools are used to simplify, analyze, and 
illuminate behaviour of a variety of chemical-reaction-based natural phe- 
nomena. Consequently, many foundational concepts of such mathematical 
models have once again achieved prominence in the mathematical litera- 
ture niEiEiinHisi- 

One such foundational concept is that of complex balancing of chemical 
reaction networks, which has been used in analysis of industrial chemical 
reaction networks [7]. In 1972, the authors F. Horn and R. Jackson showed 
that relative to every compatibility class — the invariant space where so- 
lutions lie — complex balanced chemical reaction networks necessarily have 
exactly one positive equilibrium state and that this equilibrium state is 
asymptotically stable [10]. In conjunction with M. Feinberg, they also de- 
rived necessary and sufficient conditions for a system to be complex balanced 
based solely on the reaction graph of the system [6j[8] . This work culminated 
in the Deficiency Zero Theorem and was a substantial generalization of ex- 
isting results on stability which required conservation of mass and balancing 
of forward and backwards reaction rates at equilibrium for each reaction. 

It was theorized at the time that the convergence of solutions to the pos- 
itive equilibrium state extended globally to the entire positive compatibility 
class, effectively eliminating the possibility that solutions converged to the 
boundary of the positive orthant. (Indeed, the point seemed so inextricably 
connected with asymptotic stability that in the original paper the authors 
errantly asserted that they had in fact proved just that! [10]) To date, how- 
ever, the conjecture is only known to hold for certain special cases, which 
will not be summarized here. Important work has also been done in restrict- 
ing the nature of any possible w-limit points on the boundary. In [13], the 
authors show that any w-limit point lying on the boundary is a complex 
balanced equilibrium concentration. In pQ and [3], the authors show that 
oj-limit points may only lie on certain subsets of the boundary where these 
subsets can be easily determined by the reaction graph of the mechanism. 

In this paper, we extend the stability results obtained in [5]. In that 
paper, the authors showed that for detailed balanced mechanisms with 
bounded, two-dimensional compatibility classes, solutions originating in the 
positive orthant necessarily tend to the associated positive equilibrium con- 
centration and not to the boundary. Their approach consisted of dividing 
the positive orthant into regions, called strata, and then manipulating the 
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governing differential equations of the mechanism to show that within each 
stratum trajectories were repelled from the boundary. By generalizing their 
concept of strata, we will show how their results can be extended to complex 
balanced mechanisms. 

The paper is organized as follows: in Section 2, we briefly introduce the 
relevant mathematical model for chemical reaction networks and present the 
notion of a complex balanced system; in Section 3, we extend the notion of 
strata and the linear Lyapunov functions H(x) = {a, x) introduced in [5] to 
complex balanced mechanisms, and present a few applications and examples; 
in Section 4, we give some concluding remarks including why we think this 
approach is a significant step towards proving the general Global Attractor 
Conjecture (Proposition 12. 2| h 

Throughout the paper, we will let M™ and M> denote the m-dimensional 
spaces with all coordinates strictly positive and non-negative, respectively. 

2 Background 

In this section, we outline the important concepts of chemical kinetics which 
will be needed throughout this paper. We introduce the concept of complex 
balancing first introduced in [6U511T0] and outline the relevant results of these 
papers. 

2.1 Chemical Reaction Mechanisms 

Within the mathematical literature, several distinct ways to represent chem- 
ical reaction networks have been proposed. In this paper, we will follow 
closely the complex-oriented formulation introduced by Horn et al. in |10j . 
(For examples of reaction- and species-oriented formulations, see [1] and [14j . 
respectively.) 

An elementary chemical reaction consists of a set of reactants combining 
at some fixed rate to form some set of products. We will let Aj denote the 
species or reactants of the system and define |<S| = m where S is the set 
of distinct species of the system. The set of all reactants or all products 
of a particular reaction are called complexes and will be denoted C{. They 
are linear combinations of the species and therefore can be denoted Cj = 
Y^j=i z ijAj where Zj = [zu, 2^2, • • • , Zim] G Z> . We define \C\ = n where C 
is the set of distinct complexes in the system. 

It is convenient to represent the elementary reactions of our system not 
as a list of individual reactions, but as interactions between the n distinct 
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complexes of the system. In this setting, the reaction network can be rep- 
resented as 

Ci fc -H Cj, for i,j = 1, . . . , n, (1) 

where Ci is the reactant complex, Cj is the product complex, and k(i,j) > 
is the reaction rate associated with the reaction from Ci to Cj [10]. This 
representation of a chemical kinetics mechanism will be called the reaction 
graph. 

Note that if either i = j or the mechanism does not contain a reaction 
with Ci as the reactant and Cj as the product, then k(i,j) = 0. Otherwise, 
k(i,j) > 0. The set of index pairs for which k(i,j) > will be denoted 
by 7Z and the number of such index pairs will be denoted by \7Z\ = r. 

2.2 Mass- Action Kinetics 

We are particularly interested in the evolution of the concentrations of the 
chemical species. We will let Xi = [A%[ denote the concentration of the i th 
species and denote by x = [x\ X2 ■ ■ ■ x m ] T E M> the concentration vector. 

The differential equations governing the chemical reactions system ([I]) 
under the assumption of mass action dynamics can be expressed as 

fa 

- = f(x)= J2 (■i-*)* - ' (2) 

(i,j)en 

where x Zi = YYj=i • 

Several fundamental properties of chemical kinetics systems are readily 
seen from this formulation. In particular, it is clear from ([2]) that solutions 
are not able to wander around freely in R m . Instead, they are restricted to 
stoichiometric compatibility classes, |10| . 

Definition 2.1. The stoichiometric subspace for a chemical reaction 
mechanism {7]j is the linear subspace S C M m such that 

S = span{(zj - Zj) | (i,j) G TZ} . 

The dimension of the stoichiometric subspace will be denoted by \S\ = s. 

Definition 2.2. The positive stoichiometric compatibility class con- 
taining the initial concentration xq € M™ is the set C xo = (xo + S) n M> . 

Proposition 2.1 ( [TOlfTi] ). Let x(t) be the solution to (0j with x(0) = x € 
R™ . Then x(t) G C Xo for t > 0. 
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Note that a solution x(t) of ([2|) with x(0) = xo G M> may exist only on a 
finite interval < t < T, in which case x(i) G C Xo for < t < T. Throughout 
this paper we only consider solutions to (|2|) satisfying x(0) = xo G K>0' so 
that Proposition 12.11 holds. 

2.3 Detailed and Complex Balanced Systems 

One important characteristic by which we can categorize chemical reaction 
mechanisms is the nature of the equilibrium concentrations permitted by the 
mechanism. We start by introducing two such classifications and illustrating 
how they are related. 

Definition 2.3. The concentration x* G M> is said to be a detailed bal- 
anced equilibrium concentration of ^) if 

k(i,j)(x*r = Hj,i)(x*P (3) 

for all i, j = 1, . . . ,n. A mass-action system is said to be detailed bal- 
anced for a given set of rate constants k(i,j) if every positive equilibrium 
concentration of (0|) is detailed balanced. 

Definition 2.4. The concentration x* G K> is said to be a complex bal- 
anced equilibrium concentration of |TJp if 

n n 

£fc0\0(O Bj = (^rX>fri) (4) 

for all i = 1, . . . ,n. A mass-action system is said to be complex balanced 

for a given set of rate constants k(i,j) if every positive equilibrium concen- 
tration of is complex balanced. 

Analysis of complex balanced systems is made easier by the following 
lemma. An analogous result exists for detailed balanced systems as a con- 
sequence of detailed balanced equilibria being a subset of complex balanced 
equilibria (see Theorem 3.10, |13j). 

Lemma 2.1 (Lemma 5B, |lUj). If a mass action system is complex bal- 
anced at some concentration x* G M>q> then it is complex balanced at all 
equilibrium concentrations. 

It is clear that every detailed or complex balanced equilibrium concentra- 
tion is an equilibrium concentration of ([2]) and that every detailed balanced 
equilibrium concentration is also complex balanced. It should be noted, 
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however, that not every equilibrium concentration is a detailed or complex 
balanced equilibrium concentration, as can be seen by the system 

2A 2 -A 2 Al . W 

No equilibrium permitted by mechanism (0) is either detailed or complex 
balanced. Similarly, not every complex balanced equilibrium is a detailed 
balanced equilibrium, as can be seen by 

7 \ Sfi (6) 

A 3 . 

This mechanism permits complex balanced equilibria, but no detailed bal- 
anced equilibria. 

The structure of the reaction graph is intricately connected to the con- 
ditions of detailed and complex balancing of equilibrium points; however, 
for the sake of brevity we omit such discussion here (for further details, 

see pEmnmi]). 

2.4 Known Stability Results 

In this section, we will discuss some of the known stability results for complex 
balanced systems. In particular, we state what has come to be known as 
the Global Attractor Conjecture and give several circumstances discussed in 
the literature under which it is known to hold. 
We start, however, with the main result of |10j . 



Theorem 2.1. If a mass-action system is complex balanced, then there ex- 
ists within each positive compatibility class C Xo a unique positive equilibrium 
point x* which is asymptotically stable. 

This is a very powerful result in that it gives sufficient conditions for an 
extremely predictible — and often desirable in practice — form of behaviour 
based solely on the nature of the equilibrium points. (Further relationships 
between the complex balancing condition and the reaction graph of a mech- 
anism are given in [HUE].) 

The result, however, is local in nature and insufficient to eliminate the 
possibility of an w-limit point lying on the boundary of M>q- The hypothesis 
that the unique positive equilibrium point is in fact a global attractor for 
the invariant set C Xo was first proposed by Horn in [9]. For completeness, 
we state the conjecture here as stated in [5]. 
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Proposition 2.2 (Global Attractor Conjecture). For any complex balanced 
system and any starting point xo € K>o> ^ e associated complex balanced 
equilibrium point x* of C X(J is a global attractor of C xo . 

To date no fully general proof of the conjecture exists; however, sev- 
eral restrictions on long-term behaviour of solutions and special cases under 
which the conjecture holds are known. In [13] . the authors prove the follow- 
ing restriction on the w-limit set. 

Theorem 2.2 (Theorem 3.2, [13J). For any xo £ M> °f a complex balanced 
mass-action system, the u -limit set consists either of complex balanced equi- 
librium points lying on dM™ Q or of a single positive point of complex balanced 
equilibrium. 

An important consequence of this result is that in order to show the global 
attractor conjecture holds for a complex balanced system it is sufficient to 
show that w(xo)ncM>Q = (see also Proposition 19 of [5] and the consequent 
discussion). 

Significant work has also been done recently restricting where w-limit 
points may lie on <9M> . In particular, the following concept (called a semi- 
locking set in p] and a siphon in [3]) has been used to restrict such points. 
(The set Lj is formally introduced by Definition 13.31 in Section [3.11 ) 

Definition 2.5. The nonempty index set I C S is called a semi-locking 

set if for every reaction where an element from I is in the product complex, 
an element from I is also in the reactant complex. 

Lemma 2.2 (Theorem 2.5, [1]). Consider the non-empty index set I C S. 
If there exists a xq € 1R> such that oj(xo) H Lj ^ 0, then I is a semi-locking 
set. 

Consequently, to eliminate the possibility of trajectories approaching cM> , 
it is sufficient to look at sets Lj corresponding to semi-locking sets. Since 
semi-locking sets can be determined from the reaction graph of the mecha- 
nism alone, this is a particularly useful result. 

Our work in the next section will focus on the approach used and re- 
sults obtained in [5]. In this paper, the authors use the novel approach of 
dividing the state space R> into strata, which are naturally arising parti- 
tions of M> . The authors show that within these regions, the trajectories 
of any detailed balanced system are repelled from the boundary by a linear 
Lyapunov function of the form H(x(t)) = (a,x(t)). They use this to justify 
the following partial result on global stability. 



7 



Theorem 2.3 (Theorem 23, [5J)- Consider a detailed balancing system 
whose stoichiometric substace S is two-dimensional and assume that the 
positive compatibility class C xo is bounded. Then the unique positive equilib- 
rium point x* of C Xo is a global attractor for C XQ . 

In the following section, we will generalize the methodology used in ob- 
taining this result to complex balanced systems. We feel these results repre- 
sent not only a useful tool for select special cases but a substantial theoretical 
step toward confirming Proposition 12.21 



In this section, we show how the results of [5] can be extended from de- 
tailed balanced systems to complex balanced systems. In Section 13.1} we 
introduce the necessary background material and generalize the notion of 
stratifying the state space M> - ^ n Section 1331 we derive the analogous result 
to Corollary 18 of [5] for cyclic complex balanced systems, and in Section 
13.31 we extend this result to general complex balanced systems. In Section 
I3.4l we give a few concrete results for determining global stability of complex 
balanced systems. These results are applied in Section [3.51 to a few specific 
examples. 

3.1 Permutations, Faces and Strata 

In this section, we introduce the concept of stratifying the positive stoichio- 
metric compatibility classes C X() as it has been used in the literature so far. 
We then generalize the concept in a natural way so that it can be applied 
to complex balanced chemical reaction systems. 

The idea of stratifying C X[) was first introduced in [5] with applications 
to detailed balanced systems. For such systems, the authors used the sets 



where x* is the unique positive equilibrium concentration in C Xo . The sets 
E' C TZ were chosen to contain exactly one of the index pairs (i,j) or 
out of each detailed balanced pair given in Definition 12.31 The graph of E' 
was also required to be acyclic. 

While our notion of stratification is based on that presented in [5], some 
differences arise. We consider a complete ordering of all the complexes in the 
system, rather than pairwise ordering as in ([7]), and we do not require any 



3 Global Stability 
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conditions on the reaction graph. We also keep the notion of stratification 
general by considering the state space M> rather than each C Xo . 

First of all, we will need to introduce the concept of a permutation 
operator. 

Definition 3.1. Consider the set I = {1, 2, . . . , n}. The operator u : / i— > / 
is called a permutation operator if it is bijective. Furthermore, we will 
say that the permutation operator u implies the ordering 

u(i) >- u(i + l), i = 1, . . . ,n - 1 

on the set {1, 2, . . . , n}. 

A permutation operator simply shuffles the elements of a set. To each 
such operator we can define a stratum in the following way. 

Definition 3.2. Given a permutation operator u : I i— >■ I we define the 
stratum associated with u to be 

S» = {* 6 »;'» | (^r W > for i = 1, • • • , n - 1 } (8) 

where x* is an arbitrary positive equilibrium concentration permitted by the 
system. 

This is a more general notion of strata than that given by ([7]) (i.e. for 
some systems, strata according to ([7]) are further stratified by ©); however, 
it is natural for the analysis we undertake in the remainder of this paper. 

Strata defined in this way have some nice properties, most importantly, 
that each x G M> either belongs to a unique stratum or the boundary 
separating one or more strata. 

It is also worth noting that not every permutation generates a non-empty 
stratum. For example, for a system containing the complexes C\ = 0,C2 = 
A\ , C3 = A2 , and C4 = A\ + A2 there are no points satisfying 

2 1 2 1 

since the first and last conditions imply x\ > x\ and x\ > x*, respectively. 
That is to say, for the permutation u([l, 2, 3, 4]) = [3, 4, 2, 1] we have <S M = 
(^([1,2,3,4]) = [3,4,2,1] will be our short-hand for a(l) = 3, fi{2) = 4, 
u(3) = 2, /i(4) = 1). In this paper, we will consider only those permuta- 
tion operators u which generate non-empty strata S^. (This is related to, 
although not equivalent to, the condition that E' contain no cycles in (j7J).) 
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To answer the question of global stability, we are interested in the be- 
haviour of trajectories near the boundary of the state space K> . Such 
discussion is aided by partitioning cM> into the following sets Lj. (These 
sets are defined similarly in [T], [3], and [5]. In [2], Lj is denoted Zj. In the 
standard theory of convex polytopes, Lj is referred to as the relative interior 
of a face.) 

Definition 3.3. Given an index set I C {1, 2, . . . , m}, we will define the set 
Lj to be 

Lj = {x € M> | Xi = for i £ I, and X\ > for i g 1} . 

It should be noted that according to this definition each x € cM> can be 
placed into exactly one Lj so that the Lj uniquely and completely decompose 
dR™ . 

The following result relates strata and the sets Lj. It is based on Lemma 
17 of 0. 

Lemma 3.1. // S^Pi Lj ^ then there exists an a € M. m satisfying 

on < 0, for i £ I , 

a, = 0, fori^I U 

and 

( z ^(i) - z fi(i+i),a) >0, fori = 1, ...,n- 1. (10) 

Proof. Suppose there is no a € M m satisfying ([9]) and (jlOp . By application 
of Farkas' Lemma on the index set /, this implies that there exist A.; > 0, 
i = 1, . . . , n — 1 such that 

n-l 

v = Yl A *( Z M*) ~ z Mi+l)) ( n ) 
i=l 

satisfies 

Vi > 0, for all z 6 / . , 

Vi > 0, for at least one i® € /. 

By assumption we have (1 Li ^ 0. This implies that there exists a 
sequence {x fc } C iS M such that x fc — s> x € Lj as /c — > oo. By consideration of 
the quantity (x/x*) v separately for x € Lj and the sequence {x fc } we will 
produce a contradiction. 
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Consider x G Lj. This implies Xi = for i G I. Since Vi > for is/ 
and there exists at least one iq G I such that Uj > 0, it follows that 



x \ v 



(13) 



Now consider the sequence {x fc } C converging to x. We have 

k \ 127=1 A i( z M (») _Z M(i + l)) n-l r / 1.x Z..,.n-Z..,,,,%"l A i 



x 

X" 



n-l 

n 

i=l 



X 

x" 



S ,j(i)- Z M(*+1) 



It follows from x G 5^ and Aj > that, for % = 1, . . . , n — 1, 

A, 



X 

X s 



5 ^(i)- z M(i+l) 



> 1, which implies I I > 1. 



x 



X" 



It remains to take the limit x fc — > x. The function (x/x*) v is continuous 
on M> j furthermore, it is continuous at any x G 
Xi = 0. Since v satisfies this for x fc — > x G L/, we have 



> such that f j > if 



lim 

fc-»oo V x* 



x \ v 



X" 



> 1. 



(14) 



This contradicts ()T3l) . It follows that no v satisfying (|TT|) and (|T2|) exists. 
However, the existence of such a v was a direct consequence of the non- 
existence of an a satisfying ([9]) and (llOj) . so it follows that such an a must 
exist. This proves our claim. 

□ 



3.2 Cyclic Complex Balanced Systems 

In this section, we consider the properties of cyclic complex balanced sys- 
tems. 

We start by introducing the concept of a reaction cycle as it is used 
in [10]. 

Definition 3.4. A family of complex indices {vo, V\ % . . . , v{\, I > 2, will be 
called a cycle if 

u = vi (15) 
but all other members of the family are distinct, and if 

k(uj-i,Uj) > 0, j = l,2, 
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where I is the length of the cycle. The reaction cycle associated with 
{vq, ui, . . . , v\\ is defined to be the corresponding set of elementary reactions 

C Vj _, — > C Uj , j = 1,2,..., I. (16) 

Definition 3.5. We will say that a mass-action system is cyclic if the 

system consists only of a single reaction cycle. For such a system, it will be 
understood that I = n. 

In this section, we will consider only cyclic systems. For notational 
simplicity, we reindex our single cycle of consideration to {1, 2, 1} so 
that the reaction cycle in consideration is 

Cj—i > Cj, j ' = 1, 2, . . . , I 

where Cq = C\. 

Reaction cycles are a central topic of consideration in [10] . The following 
result corresponds to equation (5-10) of that paper. (This should also be 
contrasted with Lemma 16 in [5] for detailed balanced systems.) 

Lemma 3.2. Consider a cyclic mass-action system. If the system is complex 
balanced then can be written 

t=i 

where x* is the unique positive equilibrium point guaranteed by complex bal- 
ancing and k > 0. 

Proof. Since the system is cyclic with the cycle {1,2,..., n, 1}, we can write 
© as 

dx " 

— = J>(M + l)(z m - Zi )x z > (18) 

i=l 

where i = n + 1 implies i = 1. By the assumption of complex balancing, 
from (|H) we have k(i - l,i)(x*) Zi - 1 = k(i,i + l)(x*) Zi for i = 1, . . . , n. This 
can only be true for a cyclic system if 

k(l, 2)(x*) zi = k(2, 3)(x*) Z2 = • • • = k(n, l)(x*) zn = k > 0. (19) 

Solving for each k(i, i + 1) individually, we have 

k(i, i + 1) = — — — , for i = 1, . . . , n 

v ' > (x*) Zi 

which upon substitution into (|18p yields (|17|) and we are done. □ 
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The following result allows us to rearrange the governing system of dif- 
ferential equations given by (|17p into a form which will be convenient in 
light of our conception of strata. Since we are dealing with strata, we will 
need to recall the definition of a permutation operator (Definition 13. ip . 

It is important to notice the difference between + 1) and + 1: 
the increment fj,(j + 1) is made with respect to the implied ordering given 
by the permutation //, while the increment n(j) + 1 is made with respect to 
the original ordering of the cycle. 

Theorem 3.1. Given a cyclic complex balanced system and an arbitrary 
permutation operator \i, the system J2|) can be written 



dx 
~dt 



n-1 
i=l 



V(i+i) 



(20) 



where s 



V0')+i ~ z /*Cj)- 



Proof. We notice first of all that, since the system is cyclic, we have 



Zj+1 



i=l 



i=l 
n 



i=l 



which immediately implies K ^^ s u(«) ( — ~) ' = 



i=l 



Subtracting this from (|17p . which is the form of (J2J) justified by Lemma 
we have 



, n— 1 

ax v - 

i=l 
n-1 

1=1 

n-1 



X 

X* 

n-1 



X \ z m(") 
X*/ 



^2 (( x * 



S M(j) 



3=i 



K 



E 

i=l 



X] s Mi) 



5 M(j + l) 



and the result is shown. 



□ 



It is clear from (J20j) that the vectors Y?j=i s «(i)' * = ■ ■ ■ > n ~~ ^ 
an intricate role in determining the dynamics of a system within a given 
stratum. The following result allows us to further understand the nature of 
these vectors. 
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Lemma 3.3. For every permutation operator \x and every fe = l, 2, ...,n— 1, 
there exist Xj G %<o, j = 1, 2, . . . , n — 1, suc/i i/iat 

fc n— 1 

£ s Mi) = S A J ( z Mi) ~ z Mi+i)) 
j=i i=i 

where s Mi) = z M(j)+1 - z^-) . 

Proof. Consider a permutation operator fi and fix a A; G {1, 2, ... ,n — 1}. 
Consider 

S M*0 = z M(fc)+i ~ z M(fc)- 

Clearly, there exists a t\ G {1,2, ... ,n} such that fi(k) + 1 = //(ii). We 
need to consider where n(t\) lies in the ordering implied by fi relative to 
n(k), in particular, whether (1) (J,(ti) >- /u(fc), or (2) /i(ti) -< We will 

use an iterative process on the vectors s^U), j = 1, ■ ■ ■ , k, to show that the 
case /ix(ti) >- fi(k) eventually leads us in a natural way to consideration of 
an index ti satisfying n(ti ) -< /u(fe). 

Case 1: If fJ,(t\) >~ fi(k) then s^ tl ^ is a term in the sum Ylj=i s Mj)' ^ 
follows that 

s M ( fe ) + s M(tl ) = (z M ( fc )+i - z Mfc) ) + (z Mtl)+1 - z Mtl) ) 

since + 1 = n(t±). We now repeat this process. We know that there 
exists a ^2 G {1, 2, . . . , n} such that /i(ti) + 1 = £t(*2) and, as before, either 
/i(i2) >~ n(k) or /ifo) ~< fj,(k). If fifo) ^ we add s M ( t2 ) to the cumulative 
sum (|2ip . We can continue doing this until we arrive at an index i$ for which 
KUo-l) + 1 = M&o) ^ yielding 

io— 1 

s^(fc) + s rtu) = z ku ) - z Kk) • ( 22 ) 

i=l 

We know such a terminal index exists because the cyclic nature of the system 
guarantees each index fi(ti — 1) + 1 = ^(U) >~ ^(k) is unique, so that a distinct 
vector s M ( t .) is chosen during each iteration. Since k < n and the cycle is of 
length n, this process must reach an index fj,(ti — 1) + 1 = /u(ij ) -< n(k) 
eventually. 
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Case 2: If fj,(ti ) ~< fi(k), we can interpolate ([22]) as follows: 



i — l 

s M ( fe ) + J] S M*0 = z ^(t !0 ) - z Mfc) 



( Z At(ti ) - Z /^(i l0 -l)) + ' ' ' + ( z M (fc+l) - Z M(fc)) 

~ z /*(k+i)) ( z m(*« -i) ~ Z K%))- 



(23) 



Notice that if our initial reindexing /u(fe) + 1 = ^(ii) yielded -< //(A;), 

we can take ti = t\ in the above argument. This amounts to interpolating 
s M ( fe ) = z M (*i) - z M(fc) directly. 

We return now to consideration of the entire sum Ylj=i s At(i) • Since a 
distinct vector s^ua is chosen in each application of the argument for Case 
1, we can divide this sum into those elements s^rj) considered in (|23p and 
those not. For those elements not yet considered, the same argument can be 
applied starting with the lowest remaining index, which will yield another 
sum of the form (|23p . This will remove some of the remaining vectors s^/j) 
from the sum. Since there are a finite number of complexes, this process must 
terminate at some point. Clearly, any sum of vectors of the form given in 
(I23p has non-positive integer coefficients for the terms z^(j) — z n(j+i), so that 
the existence of Xj € Z<q is guaranteed. Since \i and k E {1,2,... ,n — 1} 
were chosen arbitrarily, the result follows. 

□ 

The results to this point are sufficient to prove the following result. This 
should be contrasted with Corollary 18 of [5]. 



Lemma 3.4. Consider a cyclic complex balanced system and an arbitrary 
permutation operator fx. If fl Lj ^ then there exists an a G M< 
satisfying 

cti < 0, for i £ I 
a>i =0, for i 7 

such that {a, f (x)) < for every x £ S^. 

Proof. Since D Li ^ 0, we know by Lemma 13.11 that there exists an 
a G K<o satisfying 

cti < 0, for i E I 
cti = 0, for i 7 

such that 

( z Mi) ~ z Mi+i)' a ) - °> for i = 1, . . . ,n - 1. 
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According to Theorem 13. II we have 

/ n-l 

(a,f(x)) = / a,/c^ 

n-l 



i=l 



X WO) / X WCi+i) 



X" 



X \ z m(j) / X WU+l) 



i=l 



X" 

' i * 
a, J] s Mi) 



(24) 



For every x £ <Su, by Definition 13.21 we have 



((I 



x / x \ z ^(i+i) 

x 



> 0. 



(25) 



Now consider (a, ^}=i s «(j))- We know from Lemma [3.31 that there exist 
Xj € Z<o, j = 1, 2, . . . , n — 1, such that 

i n—1 

j=i i=i 

We also know by Lemma 1331 that (z^j) — z^+^a) > 0. Together, these 
facts imply that for every % = 1, 2, . . . , n — 1, 



/ i \ n-l 

a > s Mi) ) = J2 A i( z Mi) " z Mi+i)< a > ^ °- 

\ 3=1 / j=l 



(26) 



It follows immediately from (|24p . (|25p . (|26p and the fact that k > that 
(a, f(x)) < for every x S <S M , and we are done. 

□ 



3.3 General Complex Balanced Systems 

In this section, we extend Lemma 13.41 to general complex balanced systems. 
We follow the methodology employed by Horn et al. in generalizing from 
cyclic complex balanced systems to general complex balanced systems [10] . 

The following result extends Lemma 13.21 to general complex balanced 
systems. 

Lemma 3.5. Consider a mass- action system which is complex balanced at 
x* £ R>o- Then there exists a 5 € Z>o and K{ > 0, i = 1, 2, . . . , 5, such that 

cfoc 

— = kiXi + k 2 X 2 H h k s Xs (27) 
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where 



where the set 



..(*) 



i=i 



Z (i) 



z («) 

3 



(28) 



,(') 



, is a cycle according to Definition \3.4 



Proof. According to Lemma 6D of [10J . if a system is complex balanced at 
x* £ R™ then it can be decomposed into a finite number of cyclic subsys- 
tems which are all complex balanced at x*. This decomposition occurs with 
respect to the rate constants k(i,j), which enter linearly. This implies 
(|2l) can be written as 



f/x 
~dt 



Y 1 +Y 2 + --- + Y s 



for 5 € N+, where each Y», i = 1,2, ...,<5, corresponds to a cyclic mass- 
action system which is complex balanced at x* . Applying Lemma 13.21 to 
each of these terms yields (127j) with the difference that we must use the 
ordering of each individual cycle in transforming (jTTJ) to (|28p . □ 



We are now prepared to generalize Lemma 13.41 to general complex bal- 
anced systems. 

Theorem 3.2. Consider a complex balanced system and an arbitrary per- 
mutation operator fi. If n Lj ^ then there exists an a € M< satisying 

cti < 0, for i £ I 
ati = 0, for i £ I 



such that (a, f(x)} < for every x G S^. 

Proof. Consider a permutation operator fi satisfying S^nLj ^ 
13.11 there exists an a € M< satisfying 

cti < 0, for i 6 / 
cti = 0, for i ^ 7 



and 



for i = 1, 



1. 



By Lemma 
(29) 
(30) 



( z fi(i) ~ z ^(i+i),a) > 0, 

The form of a from (|29|) is what we need for the theorem. We now want to 
use ([30]) to determine the sign of (a, f(x)). 
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Since the system is complex balanced, by Lemma 13.51 we have 

(a,f(x)) = Ki(a,Xi) H h K<5(a,X 5 ) 

where the re, are positive constants determined by the rate constants and the 
Xj have the form (|28p . Each Xj corresponds to a cycle in the cyclic decom- 
position of the system where the i th cycle is indexed vjpi • • • , ^j. 
The overall ordering 

M(l) >- /x(2) >- ••• >- /i(n) (31) 

implies an ordering on the complex indices | , . . . , vf^ | . We can do this 
by simply removing the elements from (|31|) which do not correspond to 
indices in the set | , . . . , v^) | whilst otherwise preserving the ordering. 

Now consider a single term {a, Xj), i = 1, . . . , 5. Firstly, we reindex the 
complexes so that the relevant cycle is {l,2,...,Zj,l}. We let /Xj denote 
the permutation operator which preserves the ordering implied by \i on this 
reduced index set, after reindexing. (For example, consider a system with 
five complexes and the cycle {2, 4, 1, 2}. Consider the permutation operator 
/i([l, 2, 3, 4, 5]) = [2,5,3,1,4]. Then we reindex the cycle so that we have 
{1,2,3,1} and ^([1,2,3]) = [1,3,2] since 2 >- 1 >- 4 in the original 
ordering implied by fj,.) 

Since Xj is cyclic and complex balanced, we can apply all of the results 
used in the proof of Lemma 13.41 to get 

i=l \ 3=1 I 

where s^u^ = — z /i.,(j)- Since the ordering of the complexes corre- 

sponding to elements in the i th cycle satisfy (I3ip . we have 



X" 



x \ z Mj(j+i: 

X* / 



> 



for all x € S^. Similarly, we can apply Lemma 13.31 to show that 

k-i 

= Yl X i( z MJ) ~ z Mi+i)> a ) ^ ( 33 ) 

where Xj € Z<o for j = 1, 2, . . . , Zj — 1. This implies Ki{a, Xj) < 0. Since we 
can carry out this procedure for all i = 1, . . . , S, we have 

(a,f(x)) = Ki(a,Xi) H h ns{a,~Xs) < 
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and we are done. 



□ 



3.4 Applications 

Several global stability results follow immediately from Theorem 13.21 In 
particular, this theorem is sufficient to guarantee solutions of ([2j) do not 
approach the boundary for general complex balanced systems if they remain 
within a single stratum. This is clear because, if we take T > to be 
the final time that a trajectory x(t) enters the relevant stratum, the linear 
functional H(pc(t)) = (a, x(t)), where a satisfies ([9]), must satisfy H(x(t)) < 
H(x(T)) < for all f > T since |#(x(i)) = (a,f(x(*))) < for all t > T 
by Theorem 13.21 This contradicts the observation that, if x(i) converges to 
x* € Lj then 



If multiple strata intersect a given set Lj, however, we cannot guar- 
antee the existence of a common a satisfying (a,f(x)) < simultaneously 
within all such strata. Consequently, we cannot rule out the possibility that 
trajectories approach the boundary through creative maneouvering between 

StjI"£Ltr£L. 

This difficulty, however, does not always arise. The following results will 
complete the analysis for such systems. We have, however, purposefully kept 
the first result (Theorem [33]) general so it may be applied to systems outside 
the scope of complex balanced systems and the strata approach taken in this 
paper. We will make explicit the connection with our systems of interest in 
a later result (Corollary 13. ip . 

Throughout this section, when we say that U is a neighbourhood of K 
in M>q we mean that U is an open covering of K restricted to K> - 

Theorem 3.3. Consider a general mass-action system with bounded solu- 
tions. Suppose that for every set Lj corresponding to a semi-locking set I 
there exists an aj € M< satisfying 



and the following property: for every compact subset K of Lj, there exists a 
neighbourhood U of K in IR™ such that 



lim H(x(t)) = H(x*) = 0. 



(ai)i < 0, for i € I 
(aj)i = 0, fori&I 



(34) 



Then cj(x ) n dW 



>o — 



(a/,f(x)) < for all x G U. 
for all x G M™ . 



(35) 
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Proof. We will let |7| denote the number of elements in the set 7 (i.e. the 
number of indices i such that Xj = for x € Lj). 

We will prove that w(xo) PI cM> = by showing that cj(xo) n Lj = for 
all 7 from |7| = m to |7| = 1. This induction corresponds to the dimension of 
Lj going from (the origin) to m — 1. Since <9K> is completely partitioned 
into such sets, this is sufficient to prove the claim. 

Our inductive step will consist in showing that w(xq) Pi Lj ^ for any 
semi-locking set 7 satisfying |7| = k, 1 < k < m, implies (w(xo)PLj)\Lj ^ 0. 
This is sufficient to violate the inductive hypothesis that w(xo) PI Lj = for 
all 7 such that |7| > k. 

We take xo G M> to be arbitrary and fixed throughout the following 
induction. 

Base case: Consider |7| = m (i.e. I = S) and suppose that I = S 
is a semi-locking set. We have Lj = {0} for which K = {0} is trivially 
a compact subset. By assumption, there exists an aj £ M< such that 
(cti, f (x)} < for all x G £7 where f7 is some neighbourhood of K in 1R> - ^ 
follows that x(i) e{x£ R™ | (aj,x) < -5} for all t > and some 5 > 0. 
In other words, we can make a "cut" sufficiently close to the origin such that 
solutions do not enter the cut out area. Consequently w(xq) PI 7/ = for 
7 = S if 7 is a semi- locking set. 

Since cj(xo) PL/ = for all 7 which are not semi- locking sets by Lemma 
12.21 it follows that w(xo) P Lj = for the base case |7| = m. 

Inductive case: Consider 1 < k < m and assume that oo(xq) P Lj = 
for all |7| > k. We will prove that u(xq) P Lj = for all |7| > fc. 

Assume u(xq) Pi Lj ^ ® for some 7 such that |7| = k and 7 is a semi- 
locking set. Since every x S Lj \ Lj satisfies x £ 7/ for some 7 such that 
|7| > the inductive hypothesis implies (w(xo) P Lj) \ Lj = 0, which is 
equivalent to (w(xo) P Lj) C Lj. In order to prove the inductive step, we 
will show that assuming u;(xq) P Lj / violates (w(xo) P Lj) C 7j. 

Consider the set K = w(xo) PL/. Since trajectories are bounded by 
assumption, w(xo) is bounded, and consequently A' is a compact set. By 
the inductive hypothesis, this is a subset of Lj so that (aj, f (x)) < for all 
x € f7 where f7 is some neighbourhood of A in M> . 

Consider the linear functional 77(x) = (aj,x). By (|34|) and (|35l) . 77(x) 
satisfies: 

1. 77(x) = if and only if x £ Lj, 

2. 77(x) < for x € M^ , and 

3. j- t H{x{t)) = (a/,f(x(i))) < for allt > such that x(t) G [7. 
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Now consider an arbitrary y € K. Since y G cj(xo), U is a neighbourhood 
of y, and H (x) is continuous, we can select a sequence ( lim = oo) 

such that {x(t fc )} £ £7", lim x(t fc ) = y, and lim H(x(t k )) = H(y) = 0. 

By Property 3 of H(x(t)) given above, H(pc(t)) may not increase while 
remaining in U and, consequently, in order to approach y 6 Lj, x(t) must 
enter and exit [7 an infinite number of times. Since U is relatively open 
in R> and x(t) € M™ for all t > by Proposition 12. H we can find a 
sequence {ifc} corresponding to the entry points |x(ifc)} C (K> \ £7) (i-e. 
x(i) G C7 for i € (tk,tk))- Because trajectories are bounded and M>o \ U 
is closed, the sequence {x(tfc)} has a convergent subsequence on R>o \ U. 
We will denote this sequence {x(tfc ! )} and let y be the point such that 
limx(t/ Ci ) = y € w(xo). Since H(x(t)) may not increase for t € 

i— >oo 

and is bounded above by zero, we have > H(x(tj-)) > H(x(tk))- Since 
lim H(x(t k )) = 0, it follows that lim H(x(t k )) = H(y) = 0, so that y € L; 

by Property 1 of H(x). 

In total we have that y € o;(x )nL / n(lR>o\L T ) = Kn(R> Q \U). We recall, 
however, that U was a neighbourhood of K in K.> so that Kn(R>o\U) = 0. 
It follows that our original assumption must have been in error, so that 
w(xo) H Lj = for all semi-locking sets I satisfying |7| = k. 

Since cu(xo) n = for all I which are not semi-locking sets by Lemma 
12.21 it follows that oj(xq)DLj = if |/| = k, and our inductive step is shown. 

Since cM> can be completely partitioned into sets Lj, 1 < |/| < m, it 
follows that 



w(x ) n 



w(xq) n dK™ = 



U i« 

l<\I\<m 

and, since xq € M> was chosen arbitrarily, the result follows. 



□ 



It should be noted that, since M> decomposes completely into com- 
patibility classes, Theorem 13.31 is sufficient to guarantee persistence within 
all compatibility classes permitted by the mechanism. This stands in con- 
trast to several results in the literature which present conditions sufficient to 
guarantee persistence relative to a specified compatibility class C xo but per- 
mit other compatibility classes of the same mechanism to be non-persistent 
(see [2] and [5]). Theorem 13.31 can be easily adapted for specific compatibil- 
ity classes by considering the sets Fj = C X[) n Li throughout the argument 
rather than the sets Lj. 

We now relate Theorem 13.31 to the methodology of Section [3l 



21 



Lemma 3.6. Let Mi denote the set of permutation operators fi such that 
5^01/^0 for a fixed I. Then, for every compact subset K of Lj, there 
exists a neighbourhood U of K in ]R> such that 

uc (J s,. 

Proof. Suppose there is a compact subset K of Lj such that, for every 
neighbourhood U of K in M> , U C U^MjS^ is violated. It follows that 
there exists a sequence {x fc } C LI^MjS^ such that x fc approaches Lj as 
k — > oo. Since is compact, we may select the sequence so that x fc — >■ x 
for some specific x 6 Lj. 

Since there are finite strata, we can select a subsequence {x fc * } C for a 
fixed /i Mj; however, this implies limx fc ' = x £ S^DLi. This contradicts 

i— >oo 

M/. Consequently, our assumption was in error, and U C U^MjS^ for 
some neighbourhood C/ of in M> - The result follows. □ 

Given Lemma 13.61 and Theorem 13.31 we can see that (|35p corresponds 
to the existence of a common ai existing in all strata adjacent to a given 
set Li, which is the desired condition. In general, however, it is difficult to 
verify this condition directly. The following result gives testable conditions 
from which fl35f) follows. It also answers the question of global stability. 

Corollary 3.1. Consider a complex balanced system. Let Mi denote the set 
of permutation operators fx such that Li ^ for a fixed L. Suppose that 
for every fixed L, 1 < |/| < m, corresponding to a semi-locking set one of 
Condition 1 or Condition 2 given below is satisfied. Then the unique positive 
complex balanced equilibrium x* of C X() is a global attractor for C XQ . 

Condition 1: We will say Condition 1 is satisfied if there exists an ai & 
HL<q satisfying ft34\ ) such that, for all i = 1, 2, . . . , n — 1 and all fi 6 Mj, 



Condition 2: Consider the cycles iu\ , v\ , ■ ■ ■ , v\ >. , iq j, i = 1,2, ... ,5, 
in the cyclic decomposition of a complex balanced system according to Lemma 
\3.5[ We will reindex each cycle to {1,2,..., li, 1} and let fii, i = 1, . . . , 5, 
denote the appropriately reindexed permutation operator restricted to the 
complexes in the i th cycle. We will say Condition 2 is satisfied if there 



22 



exists an ai € M< satisfying |^| ) such that, for all i = 1,2, ... ,5 and all 
A» S Mi, 




where s w(j) = z w(j) - z w(j - )+1 . 

Proof. The proof will proceed in the following steps. We will prove firstly 
that Condition 1 or 2 is sufficient to show (a/, f (x)) < for all x 6 U^^jS^. 
We then show by Lemma 13.61 that the such systems satisfy the hypotheses 
of Theorem 13.31 so that w(xq) D cM> = 0. We then show that for complex 
balanced systems this is enough to prove the global stability of the positive 
equilibrium concentration in each positive compatibility class. 

Consider a complex balanced system. We know that w(xo) fl Li = 
for all sets L/ corresponding to non-semi- locking sets I by Lemma l2.2i We 
also know that for complex balanced systems we have w(xq) n {0} = (see 
Proposition 20 of [5], for one proof). That is to say, we need only consider 
sets Li corresponding to semi-locking sets I such that 1 < \I\ < m. 

It is clear by the proof of Theorem 13.21 that either Condition 1 (by (|33j) ) 
or Condition 2 (by (|32p ) is sufficient to prove that (a/,f(x)) < for all 
x € U ugM/^u- This implies by Lemma 13.61 that for every compact subset K 
of Li there is a neighbourhood U of K in M> such that (ai, f(x)) < for 
all x £ U. We know that solutions of ([2]) are bounded for complex balanced 
systems since, for the function 

in 

L(x) =^x i (ln(x i ) -ln(x?) - 1) (36) 

i=l 

we have ^L(x(t)) < for all i > and x G M> [TO]. It follows by Theorem 
13.31 that w(xo) fl Li = for all such sets L/. Since we have considered all 
sets Li, it follows that w(x ) n 5M™ = 0. 

Since our system is complex balanced, it follows by Theorem 12.11 that 
there is precisely one equilibrium concentration x* € M ? " in each positive 
stoichiometric compatibility class C Xo . Since there are no w-limit points on 
the boundary of the positive orthant, by Theorem 12.21 it follows that the 
only w-limit point is the positive equilibrium concentration. It follows that 
x* is a global attractor for C XQ and we are done. □ 

Since Condition 1 implies Condition 2 by Lemma 13.31 but the converse 
does not necessarily hold, it is typically preferable to check Condition 2. 
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In the following section, our approach will be to define a set P of vectors 
Ylj=i s /i,(i)' * = 1) • • • > ^> k = • • • ' k ~ 1) an d check Condition 2 relative to 
this set. 

The following result corresponds to Corollary 4.5 of [2]. It is a general- 
ization of Theorem 23 of [5] (stated Theorem 12.31 here) to complex balanced 
systems. 

Corollary 3.2. Consider a complex balanced mass-action system whose sto- 
ichiometric subspace S is two-dimensional. Then the unique positive com- 
plex balanced equilibrium x* of C xo is a global attractor for C xo . 

Proof. With application of Corollary 13.11 the proof follows identically to the 
proof of Theorem 23 contained in [3]. We also notice that since trajectories 
of any complex balanced system are bounded by L(x(t)) < L(xo) for all 
t > 0, we may remove the assumption of boundedness. □ 

3.5 Examples 

In this section, we present two examples. The first is an example where 
Corollary 13. II can be applied while the second is an example where it cannot. 

Example 1: The following example is given in [2] as an example of a 
three-dimensional complex balanced system for which a general method of 
guaranteeing global stability is not known. The system considered is 

Ai A 2 Ai + A 2 A1 + A3. (37) 

We assign C\ = Ai, C2 = A2, C3 = Ai + A2, and C4 = Ai + A3, and 
X\ = [Ai], X2 = [A 2 ], and £3 = [.A3]. The system is complex balanced at 
all equilibrium concentrations so we need not consider the rate constants. 
The compatibility class C xo = M> is three-dimensional and the only non- 
trivial semi-locking set is I = {1,2} so that we need only consider the set 
Lj corresponding to this index set. 

We will show that all strata such that n £{1,2} 7^ have a common 
ai G M™q satisfying ()34|) and Condition 2 of Corollary 13.11 There are six [i 

such that n £{1,2} 7^ $ : 

(1) M ([l, 2, 3, 4]) = [2, 4, 1, 3] (4) /x([l, 2, 3, 4]) = [2, 1, 4, 3] 

(2) fi([l, 2, 3, 4]) = [4, 2, 1, 3] (5) M ([l, 2, 3, 4]) = [1, 2, 4, 3] 

(3) n([l, 2, 3, 4]) = [4, 1, 2, 3] (6) M ([l, 2, 3, 4]) = [1,4, 2, 3]. 

Since the vectors Y^j=i s ^i{j) are *he vector coefficients of the bracketed 
strata terms in (I20p . it is instructive to rewrite the system of differential 



24 



equations ([2]) implied by the network (|37p according to Theorem 13. 11 (This 
analysis is not, however, required to determine the set of all admissible 
vectors Yli=i s n,(i)-) We will carry out the analysis for one stratum and 



leave the rest as an exercise. 



The first stratum is given by 



x € 



»3 



X 2 Xl X 3 X\ X\ x 2 



(38) 



Since the system can be decomposed into the cycles {1,2,1}, {2,3,2}, and 
{3,4,3}, according to Lemma [331 and Theorem 13.11 the system (J2j) can be 
written 



dx 
~dt 



1 



X 2 Xi 



+ «3 



+ K 2 



1 

-1 



1 






x 2 Xl x 2 



Xi x 3 Xl 



J, ■ > 



X 



X2 



(39) 



where x* 



[xl, x 2 , x^] T is the unique positive complex balanced equilibrium 
point and k\, k 2 and K3 are positive constants determined by the rate con- 
stants. In Sfj, the bracketed terms of (I39p are strictly positive so that the 
sign of (ai, f (x)) is determined by the vector terms alone. Consider a vector 
ai € M< satisfying (i34"|) for which 

ai = Ai(-1,0,0)+A 2 (-1,-1,0), Ai > 0,A 2 > 0. 

For any such 07 we have (a/, f(x)) < 0, which is sufficient to show the linear 
function H(x(t)) = (aj,x(i)) repels trajectories from the set -^{1,2} m the 
first stratum. 

A similar analysis can be carried out in the five other strata. Removing 
repetition, the set of admissible vectors Ylj=i s n(j) ^ s 



P 



1 




-1 


-1 




1 











1 








1 



1 



Since ai = (—1,-1,0) satisfies (aj,v) < for all v £ P, we have that 
(a/,f(x)) < for all x G where 5 M is such that n £{1,2} It 
follows by Corollary 13.11 that x* is a global attractor for C 



X() 



^>0 



and we 
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are done. 



Example 2: Consider the system 

Ai ^ 2A 2 
t I 
A 2 + A 3 <- Ai+A 2 . 

We assign C\ = Ai, C% = 2A2, C3 = At + A2, and C4 = A2 + A3, and 
x\ = [Ai], X2 = [A2], and X3 = [A3]. The system is complex balanced at 
all equilibrium concentrations so we need not consider the rate constants. 
As in the previous example, the compatibility class C XQ = K> is three- 
dimensional and the only non-trivial semi-locking set is I = {1,2} so that 
we need only consider the set Lj corresponding to this index set. 

We will show that there is no common aj satisfying (|34p and Condition 
2 of Corollary 13.11 for all strata such that 5^ n £{1,2} 7^ 0- Notice that 
since Condition 1 of Corollary 13.11 implies Condition 2 by Lemma 13.31 this 
is sufficient to show that neither condition is satisfied. There are five \i such 
that n L {1 2} + 0: 



(1) ^([1,2,3,4]) = [1,4,2,3] 

(2) ^([1,2,3,4]) = [4,1,2,3] 

(3) ^([1,2,3,4]) = [4,2,1,3] 



(4) /x([l,2,3,4]) = [1,4,3,2] 

(5) /x([l,2,3,4]) = [4,1,3,2]. 



Again, we carry out the analysis for the first stratum and leave the rest 
as an exercise. The stratum is given by 

2 



x G 



»3 
So 



El > El . El > 

x l x 2 x 3 



•Z'2 



> 



Xl X2 



Xl 



(40) 



\X 1 « X o 1 X 



*1T ; 



1; A 2i "°3J 



is the unique positive complex balanced equilibrium 



where x* 
point. 

Since the system can be decomposed into the cycles {1,2, 1} and 
{1, 2, 3, 4, 1}, according to Lemma [331 and Theorem 13. 1\ the system ([2]) can 
be written 



(ix 

M =K1 



1 



2 





+ 



X2 X3 



■"»"■ 




-1 

2 




( Xl 







\x* 


1 







« 


-1 





^2 X 3 



3, ■■■■ 



El 



Xl x 2 



X] 



V 



(41) 



2(3 



where k\ and K2 are positive constants determined by the rate constants. 
Consider a vector ai G M< satisfying ()34|) for which 

a 7 = Ai(-2,-l,0) + A 2 (0,-1,0), A x > 0,A 2 > 0. 

For any such otj we have (q/, f (x)) < for all x € S^. 

A similar analysis can be carried out in the four other strata. Removing 
repetition, the set of admissible vectors 2~^ =1 s fj,(j) is 
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1 




1 




1 
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-1 


2 




1 









-1 




-2 




-2 




1 







-1 




-1 




-1 









-1 







there 


is no ai 


satisfying 


(1341) such that (a/,v 


) < for all v 



and consequently Corollary 13.11 cannot be applied. 



4 Conclusions 

In this paper, we have presented several results which extend the breadth 
of globally stable complex balanced systems (Theorem I3.3[) and contributed 
several important theoretical steps toward proving the Global Attractor Con- 
jecture (Proposition 12. 2\) . In particular, we have extended the notion of 
stratification of the state space M> to a form applicable to complex bal- 
anced systems (Definition 13. 2\i and shown that, while trajectories lie within 
a particular stratum 5^, they are necessarily repelled from any adjacent 
set Lj, the relative interior of a face, by the linear Lyapunov function 
if(x(t)) = (a,x(i)) (Theorem E2J). 

The main result of this paper, Theorem l3.2( is limited in that it prevents 
trajectories from approaching the boundary from within a single stratum 
but does not necessarily prevent convergence to the boundary for trajectories 
which continually jump between strata. This is because the linear Lyapunov 
function H(x(t)) = (a, x(i)) which repels trajectories from the boundary is 
specific to each stratum S^. Even in cases where every H(x(t)) = (a,x(t)) 
individually bounds trajectories away from a common set Lj, it cannot 
be ruled out without further work that trajectories approach Lj through 
creative maneouvring between strata. We presented one result (Corollary 
13. ip where this difficulty could be resolved and therefore global stability 
could be shown. 

It does not seem probable, however, that a trajectory should be permit- 
ted to approach a set Lj even without a common a satisfying ([9]) and (fT0|) 
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for all such that S^DLj ^ 0. This is especially true in light of the restric- 
tions placed on trajectories by the global Lyapunov function L(x(t)) given 
by (|36[) . As such, we feel the solution to the global attractor conjecture lies 
in a fuller understanding of the geometry of the strata 5„, the relationship 
to the global Lyapunov function L(x(t)), and how the functions H{x(t)) 
relate to one another for different strata. 

Acknowledgements: We would like to thank the reviewers for their 
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